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Abstract: 

A calculation of the general relativistic apsidal or perihelion precession of planet Mercury is presented, based on the Schwarz schild 
gravitational metric field tensor and the mathematical perturbation theory and yielding the same result as in Albert Einstein's 1916 
paper. Both the classical and the general relativistic contributions to Mercury's precession result in breaking Noether's symmetry of 
Kepler's first law by which a planetary closed orbit exactly repeats itself periodically, thus the Laplace-Runge-Lenz vector is not 
anymore conserved. A clear calculation of the general-relativity corrections to both Newton's gravitational potential and force is also 
shown. 
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I THE SCHWARZSCHILD METRIC 
AND THE GEODESIC EQUATIONS 

An exact solution to the 6 independent Einstein field equations 
describing an elliptical orbital motion is the Schwarzschild 
metric field tensor, whose squared line element is 

( 1 ) 

ds 2 = c 2 (l-^jdt 2 -(l- 2 -±y dr 2 - r 2 de 2 - r 2 sin 2 6d4> 2 , 
where 

fi = GM/c 2 . W 


li = 0 , 1 , 2 , 3 . 

By means of (5) and (6),we get the 4 spacetime geodesic 
equations [6] 



6+ -rO — sin0cos0</> 2 = 0, (9) 

r 


r 2 sin 2 0(f) = h. 


Defining the squared Lagrangian for massive bodies as 

^ = guv** 1 * 1 ' (3) 

= c 2 . 

where 

(4) 

X ^ = dx^/dT, 

and r is the affine time parameter, we get, using (1), 


(5) 




^6 2 + sin 2 6 

The Euler-Lagrange equations are 


d ( 

' dL\ dL 
dxV ) ~ dxV ~ 

(6) 

dT \ 



where h, the angular momentum per unit rest mass, and /c, the 
total energy per unit rest energy, are constants. 

For a celestial massive body moving in the plane 


6 = 7t/2, 

the set of the 4 spacetime geodesic equations reduces to 



where 


r 2 — r 2 <j> 2 = c 2 . 


(14) 
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Substituting (11) and (13) into (14), we obtain 


, h 2 / 2 GM\ 2 GM - , N 
r 2 + ~{ 1 — )- =c 2 (* 2 -l). 


(15) 


c z r / r 

In (15), letting 

r —► oo and r = 0 , 

one gets 

k 2 = 1 . 

By substituting 

Jr Jr J</> /i Jr 
Jr d(f) Jt r 2 d(j) 

into (15), we get 


(17) 


( 18 ) 


(--V 

\r 2 d(f)) 


h dr\ 2 h 2 x 2GM IGMh 2 

+ - = c 2 (k 2 - 1 ) +-+ ■ 

r z r c z r J 


r»2 j-3 


With the substitution 

m = 1/r, (19) 

(18) becomes 

.2 


(du \ 2 C 2 2 

U) +m - 1} 


( 20 ) 

i 2GMu i 2 GMm 3 
+ _ P _+ c 2 ' 


By differentiating (20) with respect to (f) , we obtain 

a differential equation describing the elliptical motion of a 

massive body in Schwarzschild geometry and given by 

d 2 u GM 3GM 2 

^ + “=1F + —- < 2 » 

For a circular motion, 

r = constant, 
r = r = 0 , 


thus, 


u — 1 /r = constant, 


and ( 21 ) becomes 

GM 3GM 2 ( 22 ) 

“=-r+—«■ 

II EQUATION OF MOTION FOR PFANETARY ORBITS 
IN NEWTONIAN GRAVITY 

The classical-mechanics Lagrangian of a celestial object of mass 
m influenced by the central gravitational field force of a massive 
body M in the equatorial plane 
6 = 7t/2 


L = im(f 2 + r 2 ^ 2 ) + ^^. ( 23 ) 


From the Euler-Lagrange equations, we get 

r 2 <t> = h, (24) 

.. h 2 GM 

r= 7*~ — 

Additionally, the motionally constant total energy is given by 


(25) 


E = ^m(r 2 + r 2 (p 2 ) — • ( 2 


By plugging (19) and (24) into (26) and using 
du _du d(p _ h du ^ 7 ) 

dt d(p dt r 2 d(f) 

we obtain the Newtonian equation of motion for general 
elliptical planetary orbits 


d 2 u 


GM 


d(j) 2 h 2 


(28) 


III SOEUTION TO THE NEWTONIAN ORBITAL 
EQUATION 

The function 

<»> 

is a particular solution of the nonhomogeneous linear second- 
order differential equation (28), identified as the reciprocal of 
the semi-latus rectum , and represents a circular orbit. 

The general solution of the homogeneous equation of (28) is 
u B =Ccos(<f>-<!>„), (30) 

where C is a constant of motion. 

The general solution of (28) is 

W = W A+ U B 


= + Ccos(<£ — <t>„). 


(31) 


Defining the eccentricity of the orbital trajectory as 
„ Ch 2 


GM 


(32) 


and knowing that the eccentricity is defined in function of the 
magnitude of Laplace-Runge-Lenz vector A as 


where 


we obtain 


mk 


k = GMm, 


(33) 


(34) 


C= ,A 2 • (35) 

n m 


is given by 


Setting 
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4>« = 0. 


one gets 


u = ^^-(1 +ecos(f )). (36) 


h 2 


IV SOLUTION TO THE EINSTEINIAN ORBITAL 
EQUATION 

In (21), 

3 GMu 2 ^ GM ( 37 ) 

<<[ W 

In fact 

3 GMu 2 GM _ ^GM 2 GM 

c 2 f ' 2 c 2 r ( rv ) 2 (38) 

= 3 v 2 /c 2 , 

which is in the case of Mercury, the fastest solar planet, 


Let 


1(T 7 . 


3 (GM) 2 

€ = ' «L 

h 2 c 2 


(39) 


Thus, (21) can be written as 


GM ( hW 
U +U =* +e {™ 


(40) 


(41) 


(38) is a non-linear and nonhomogeneous second-order ordinary 
differential equation where 

€ = « 1 . 

Therefore, (40) can be solved by means of the well-established 
perturbative method, and the general solution of (38) has the 
form 

U = e°u 0 + eti! + 0(e 2 ). 

= u 0 + eu x + 0(e 2 ). 

By setting e equal to zero in both (40) and (41), we get the 
Newtonian classical orbital equation of motion whose general 
solution is (36),i.e., 

GM 

u = u 0 = —-(1 +ecos</>). (42) 

hr 

Consequently, as a first-order approximation, the Newtonian 
solution u 0 is the zeroth-order solution of the general- 
relativistic orbital equation of motion ( 21 ). 

By differentiating (41) twice and substituting the result into 

(40), we get (43) 

„ OM h 2 u 2 2x 

u 0 + Uo - -r-z 4 - e{u x 4 - u\ -) 4 - 0(e 2 ) = 0 

h GM 


€° ( u” 4- Uo — )f e(u" + U\ — ^-^) + 0(e 2 ) = 0 • ( 44 ) 

h GM 

Equating both sides of (44), we get 

„ GM _ (45) 

u" + U°- ^ =0 ' 

u" + Ui - = 0 - (46) 

GM 

By substituting (42) into (46), one gets 

u" 4- u\ = (1 + 2 ecos(j> 4- e 2 cos 2 <f>) 


or 


(48) 


n m I 2 \ 2me me 2 

u x +u i = p(l + 2 e ) + —fa2~ cos< P + ~2h2 cos ^ ’ 

which can be written as 

u i + u i = D 4- E cos<f> -}- F cos2<t >. (49) 

By choosing the general solution of (49) to be of the form 
u\ = A 4- Bcbsind) 4- Ccos2d) 

(50) 

and differentiating u x twice, we obtain 

u[ 4- u\ — (A) + (2 B)cos(f) 4- (-3C)cos2(j). (51) 

by simple comparison between (48) and (51), one gets 

. GM / i . 

A =j ? {l + -e 2 ), (52) 


D = 

C = — 


GMe 


Therefore, 


W ’ (53) 

GMe 2 

* ( 54 ) 


GMe 2 


(55) 


and 


Ul ~ + <t>sin<j> — 2 cos2<p , 

the general solution up to the first order in e of (40) is 
u = u Q + G ^ e [1 4- e(t>sin<t) 4- e 2 (^ - ^cos2(j))}. (56) 

By substituting (42) into (56), we get 

GM GM GMe 

«- ir+ir «~*+ p- 

GMe e( p slno _|_ GMe e 
+ h* h? 2 (57) 


GMe 7 . 


6/i 2 


e cos2<f> • 


All terms with e can be neglected except 


G K e e<t>sin<t>, (58) 


since the value of <p continuously increases with time as the 
planet moves on its elliptical path. Hence, (57) becomes 
u = (1 4- ecoscj) 4- ee<t>sin<j >). ( 59 ) 


or 
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Using 

cos[</>(1 -6)] = cos</>cos 6<£ + sin sin6</> 

% cos </> + 6(f) sin 4> for 6 1, 

(56) becomes 

GM . ,,,, ( 61 > 

« = ^ il+eco»[^(l -e)]]. 

V APSIDAL ANGLE AND CLOSED ORBITS 

In (25), 

GM 

~~r r ’ 

represents an attractive gravitational central force per unit mass. 
For any per-unit-mass conservative central force, (25) is 
generalized to 

?-£ = ft(r). 


Let 


( 63 ) 


r(t) = a + p(t), 
where a is the radius of a circular orbit and p is a small 

perturbation such as 

p(t )« a. (64) 

Substituting (63) into (62), we get, up to the first-order 
perturbation, 

h 2 /„ „p\ ^ dF r (a) 


For a circular orbit, (62) becomes 

h 2 „ . v 

- ^3 = F r (r). 

Plugging (66) into (65), we obtain 
-SFAa) 


( 65 ) 


dr 


( 66 ) 


P + 


F r (a) 


P = 0, 


where 

F r (a) = . ( 67 ) 

dr 

For a simple harmonic motion, the term in square brackets is 
positive , i.e., 

F r (a) + — Fr(o) < 0, ^ 

3 

and stands for the squared orbital or angular frequency. Thus, 
the circular orbit of radius a is stable to small perturbations. 
Additionally, the nearly circular orbit of the radial distance r(t) 
is elliptically bounded since r(t) oscillates between a fixed 
maximum value and a fixed minimum value called apsides, 


where, in the case of the Sun, the minimum apsis is the 
perihelion and the maximum apsis is the aphelion. 

The orbital time period of oscillation is given by 
r -fi 


T =2tt 


-3 Fr(a) 


- FAa) 


a 


( 69 ) 


An attractive central power-law force has the form 
F r (r) = -Kr n , ( 70 ) 

where n is an integer and 


K > 0 . 

Substituting (70) into (68), we obtain 

-*«!«!»“”< 0 . ( 71 ) 

3 

(71) yields 

n>- 3 . (72) 

The apsidal angle,cp, is the angle between two consecutive 
apsides, a periastron and an apastron, thus the half-cycle angular 
half-period. Additionally, the time taken by the orbit to sweep 
one apsidal angle is clearly half of the orbital time period. 

Since the angular velocity 60 or 0 is the angle through which the 
massive object turns in one second, the apsidal angle is given by 

<p- — u. ( 73 ) 

2 


But the angular momentum per unit rest mass, h, is a constant 
of motion , and the radial distance r is nearly constant for nearly 
circular orbits; hence we can write, using (66) and (24), 



Fr(a) 

a 


( 74 ) 


By means of (69) and (74), the apsidal angle per half-revolution 
becomes 


<p = n 


aF r (a) 
F r (a) 


( 75 ) 


or, 


<P = 


71 


( 76 ) 


(3 + n)V 2 

A closed orbit is an elliptical orbit that exactly repeats itself 
within every orbital time period, and whose apsidal angle is a 
rational fraction of 271 radians, which is the case of an inverse- 
square force law , where n = — 2 , such as Newton's 
gravitational law and the electrostatic law, and a linear force 
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law, where n = 1, such as Hook's law. Hence, (72) and (76) are 
in total agreement with Bertrand's theorem where only 
elliptically bounded orbits influenced by the inverse-square law 
and the radial-harmonic-oscillator law are stably closed. The 
converse is also true since all closed orbits are elliptical. 
Therefore, the apsidal angle for an elliptical orbit obeying an 
inverse-square central-force law is tt radians. 

VI GENERAL-RELATIVISTIC PRECESSION OF 
MERCURY 

Apsidal or perihelion precession of a planetary orbit is the 
rotation of the direction joining the apsides, which is the 
elliptical major axis, around an elliptical focus, which is the Sun 
in the case of the solar system, in the plane of the orbit. 

Since Kepler's first inverse-squared force law yields closed 
orbits where the perihelion or aphelion remains constant in time, 
it fails to indicate any planetary precession, which is certainly in 
contradiction with the observational fact that the perihelion, the 
aphelion, and the whole elliptical orbit rotationally precess in 
the direction of the planet's rotation round the focal Sun. 

For Mercury, the observed dynamical precession is found to be 
575"per century, where 532"per century are nonrelativistic 
and due to the gravitational disturbances of the Sun as well as 
6 outer planets—Venus, Earth, Mars, Jupiter, and Saturn—and 
can be estimated in classical mechanics. But the remaining 
43"per century cannot be account for in Newton's theory of 
gravity [28] [4]. 

In 1915 and 1916, Albert Einstein, by means of his theory of 
general relativity, could accurately explain and solve the 
centurial 43"discrepancy, which is found to be purely general 
relativistic and due the elliptical motion of Mercury in the 
gravitationally curved spacetime of the Sun. 

In fact, the general-relativistic orbital solution (61) yields an 
angular radian period 

27T (77) 

1 — e 1 

which is 2 times the apsidal angle with 

Ojr 

«27r(l + e), (78) 

thus an irrational fraction of 2 tt, yielding an open elliptical orbit. 


Since 2n is the value of the apsidal angle for a complete orbital 
cycle of a Keplerian closed orbit, the additional term 2 ti£ 
represents the general relativistic correction to the Keplerian 
apsidal angle for a full revolution. Hence, by using (39), the 
general relativistic contribution to the planetary precession per 
revolution is 


Acp = 2ne = 


6? t(GM) 2 
h 2 c 2 


( 79 ) 


For nearly circular orbits, such as Mercury's orbit, the radial 
distance r is almost constant and equal to the semi-latus rectum, 
thus (25) becomes 

h 2 = GMr (80) 

= GM( 1 -e 2 )a. 

From Kepler's third law, 

2 4;r 2 a 3 (81) 

T 2 =- 

GM 


Plugging both (80) and (81) into (79), one gets [1] 


Acp = 247r 3 


T 2 c 2 (l-e 2 ) 


( 82 ) 


which is exactly the same as in Albert Einstein's 1916 paper 
titled "The Foundation of the Generalised Theory of Relativity." 
For Mercury, the orbital period in Earth days is 


T = 88 days, 
the semi-major axis is 

a = 5.8 x 10 10 m, 
and the eccentricity is 

e = 0.2, 

and using 


( 83 ) 


( 84 ) 


( 85 ) 


Mq = 2 x 10 30 kg (86 ^ 

as the solar mass, the general relativistic precessional 
contribution per revolution amounts to 0.104". 

In an Earth century, Mercury makes 415 revolutions, yielding a 
centurial general relativistic precession equal to 43". 

In celestial mechanics, every Keplerian orbit is associated with 
two vectors that are constant throughout the planetary motion, 
thus conserved: the orbital- angular-momentum and the Laplace- 

Runge-Lenz vector, which is defined by 

-» f 

A = p x L — mk - 
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where e is the eccentricity vector. 

The LRL vector always points towards the perihelion along the 
direction of the symmetry axis of the orbit, lies in the plane of 
orbital motion, and, in consequence, fixes geometrically the 
direction of the major axis—which passes through both the 
perihelion and the aphelion—and yields a Noetherian symmetry 
by which a Keplerian orbit repeats itself after every orbital 
period of time. 

Any departure from the inverse-squared force law breaks up the 
symmetry and the conservation of the Laplace-Runge-Lenz 
vector; as a result, the LRL vector starts to rotate slowly round 
the Sun, causing the elliptical orbit of planet Mercury to precess. 


VII CORRECTIONS TO NEWTON’S 

GRAVITATIONAL FORCE AND POTENTIAL 

The nonrelativistic or classical contribution to the precession of 
a planetary orbit follows Newton's inverse-square force law, 
thus does not result in any correction to the gravitational force, 
unlike the general relativistic precessional contribution, which is 
given by (79) or (82) and, when purely considered, leads ,by 
means of (78),(79),and (80), to a purely relativistic precessional 
apsidal angle per cycle written as 
29 = 27 t (1 + e) 


= 2tt(1 + ^), 


( 88 ) 


and to a correction to Newton's gravitational force per unit mass 
given by 


GM± 


( 89 ) 


such that the classical gravitational force per unit mass in the 
case of the Sun-Mercury system becomes [12] 


CM e _ GM q 1 


( 90 ) 


GM q 3 ( GMq ) 2 


_ GM o - 3 GM q J^_ , 


r 4 c 2 


which can be demonstrated using (75). 

In fact, the derivation of (90) with respect to r yields 


Fr = 2 cm o-i-12 GM q ft 2 . 

r 3 r*c 2 

Also, with a simple calculation, we get 


T - 2 (‘ + ^)( ,+ ^) 


6/i 2 


3 /i 2 ' 


Since 


— « 1 , 


(92) becomes 

rFr o/l , 6/l 2 


7^ =" 2 ( 1+ ^)( 1 -^) 


Taking into account that 

18 /j 4 


rV 


« 1. 


( 91 ) 


( 92 ) 


( 93 ) 


( 94 ) 


( 95 ) 


(94) can be written as 


'A=- 2 (i+m (96) 

p v r 2 <r 7 


Plugging (96) into (75) for nearly circular planetary orbits, we 
get the general relativistic apsidal angle for a half-revolution 
<P = 7r(l + l^-). ( 97 ) 

multiplying (97) by 2, we obtain the general relativistic apsidal 
angle for a complete revolution 

2<P = 2tt(1 + 4^-) 

v r £ <r ' 




( 98 ) 




which is identical to (88) with a full-revolution general 

relativistic precession equal to 

67TGM (99) 

* 

rcr 

which is certainly also identical to (76). 

Defining the classical gravitational potential to be the work done 
by a conservative gravitational-field force to bring a unit mass 
from infinity, where the gravitational potential approaches zero, 
to a point at a radial distance r, we can write 

r 




( 100 ) 


and substituting (80) and (90) into (100), we obtain 
GMq 3 (GMq) 2 (101) 


V= - 




where 
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3 (GMq) 2 
2 rV 


( 102 ) 


is the general-relativity correction to Newton's gravitational 
potential. 

Unlike Newton' theory of gravity, Einstein's general relativity 
accurately estimates not only the additional precession due to 
curved spacetime but also the classical precessional contribution 
caused by other planets, which is calculated within the frame of 
general relativity in the case of weak gravitational fields where 
the curvature of spacetime is of very small value, yielding 
Newtonian gravity. 
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